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Abstract. By using the lattice Monte-Carlo simulation, we investigate the finite temperature 
chiral phase transition at color SU(3) gauge theories with various species of fundamental 
fermions, and discuss the signal of the (pre-)conformality at large Nf (num. of flavors.) via 
their comparisons. With increasing Nf, we confirm stronger fermion screenings resulting from 
a larger fermion multiplicity. We investigate a finite T step-scaling which is attributed to 
the uniqueness of the critical temperature (Tc) at each Nf, then the vanishing step-scaling 
JL ■ signals the emergence of the conformality around Nf ~ 10 — 12 . Further, motivated by the 

recent functional renormalization group analyses, we examine the Nf dependence of Tc, whose 
, vanishing behavior indicates the conformality at NJ- ~ 9 — 10. 



1. Introduction 

Conformal invariance is anticipated to emerge in asymptotically free non-Abelian gauge theories 
when the fermion species (flavor A'^^) exceeds a critical value Nf = Nf* [Il|2]. The approach to 
the conformality from below is in principle associated with a pre-conformal (walking) behavior of 
the running coupling, which has been advocated as a basis for strongly interacting mechanisms 
of electroweak symmetry breaking 

Recent lattice studies [1] focused on the computation of Nf* and the analysis of the conformal 
window itself, either with fundamental fermions or other representations. Among the many 
interesting results with fundamental fermions, we single out the observation that the color SU(3) 
gauge theory with Nf = 8 is still in the hadronic phase [5l|6], while Nf = 12 seems to be close 
to the critical number of flavors, with some groups favoring conformality [SHH], and others chiral 
symmetry breaking [10). 

In order to attack the walking and the conformal dynamics, it is more informative beyond a 
fixed Nf to investigate the vanishing or reducing chiral dynamics with increasing varying Nf. To 
this end, we investigate the Nf dependences of the chiral phase transition at finite temperature 
(T) based on our recent works [5l|Tl]: The vanishing (reducing) finite T step-scaling which 
attributes to the uniqueness of Tc at each Nf signals the emergence of the (pre-)conformality. 
Further, motivated by the recent functional renormalization group (FRG) studies |12) . we 
examine the Nf dependence of whose vanishing (decreasing) behavior indicates the (pre- 
) conformality. This thermodynamic lattice study for the large Nf non-Abelian gauge theory has 
played a crucial role to extract a notion of more strongly interacting non-Abelian plasma |13j . 
and it is expected to provide a new connection between the lattice and the Gauge/Gravity 
duality 



2. Simulation setups 

Simulations have been performed by utilizing the publicly available MILC code [15]. We use 
an improved version of the staggered action, the Asqtad action, with a one-loop Symanzik [16] 
and tadpole [17j improved gauge action. The tadpole factor uq is determined by performing 
zero temperature simulations on the 12^ lattice, and used as an input for finite temperature 
simulations. To generate configurations with mass degenerate dynamical flavors, we have used 
the rational hybrid Monte Carlo algorithm (RHMC) [18j . 

Our observables are the chiral condensate and the Polykov loop 

N 1 / ^* \ 

where Ng (Nt) represents the number of lattice sites in the spatial (temporal) direction and 
f^4,fx is the temporal link variable. The other important observable is the ratio of a scalar and 
a pseudo-scalar susceptibility [E] , 

^Xo_ ^ d{'4)'4))/dm ^ Xconn + Xdisc ^2) 

where a^^conn = -Nf{TT[{MM)-^])/{AN^,Nt) and a^Xdisc = N}[{Ti[M-^f) - 
(TiM^^)"^]/ {IQN^Nt). Here, Rj^ ~ C'(l) indicates the scalar and pseudo-scalar degeneracy 
attributed to the approximate chiral restoration, while R-,^ ^ 1 indicates the chiral symmetry 
breaking [5t[T9]. 

3. Results 

We evaluate the thermalized ensemble averages of the chiral condensate (PBP) and Polyakov 
loop (PLOOP) for various lattice couplings /3l, lattice sizes with the finite T set up A''^ ^ Nt, 
and the number of flavors Nj . All results have been obtained by using a single value for a lattice 
bare fermion mass am = 0.02. Then we locate the lattice bare coupling /3j^ associated with 
the chiral crossover which is signaled by the drastic decrease (increase) of PBP (PLOOP) as 
a function of In practice, the ratio of the scalar and pseudo-scalar susceptibility Rt^ gives 
a stronger signal of the chiral crossover, owing to its renormalization invariant property. In 
table [U we summarize the obtained critical lattice couplings as a function of {Nf,Nt). We have 
confirmed the (approximate) asymptotic scaling for the normalized critical temperature Tc/Al 
varying Nt at each Nf, where Al is so-called lattice Lambda. This indicates that our (3^^ have 
been determined near to the continuum limit [20j- 

Table 1. Summary of the (pseudo) critical lattice couplings 13^^ for the theories with 
Nf = 0, 4, 6, 8, am = 0.02 and varying Nt = 4, 6, 8, 12. The result has partially been 
extracted from our recent paper [5l[TT]. All results are obtained by using the action with the 
same level improvements. 
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7.35 ±0.1 


7.88 ± 0.05 


8.20 ±0.1 
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5.60 ±0.1 


5.89 ±0.05 


6.10 ±0.1 




6 


4.65 ± 0.05 


5.05 ±0.05 


5.2 ±0.05 


5.55 ±0.1 


8 




4.1125 ±0.0125 


4.20 ±0.1 


4.34 ± 0.04 



We shall now discuss the emergence of the (pre-)conformality by using our critical lattice 
9l ~ V^O/^L collection. The uniqueness of a critical temperature at each Nf, = 



Nt a(/5L^) = Nl a(/5L^') with Nt / Nl gives a thermal step-scaling Af3j^^ = /3l' - /3l A 
vanishing (decreasing) A/]-^ can be the signal of the (pre-)conformality. 

Our thermal step-scaling is a function of A'^^ and two lattice temporal extensions, 
APj^^Nf, Nt, Nl), and tends to be smaller with increasing Nf. We here estimate the number of 
flavor satisfying A/3^{Nj) = by extrapolating our f3j^ collection into the larger flavor region. 

To this end, we plot the (pseudo) critical lattice coupling gj^ = y/lQ/ f3j^ as a function of Nj in 
Fig. [H which gives an extension of Miransky-Yamawaki diagram [2T1 to finite T cases. 

Let us first pick up the lattice critical couplings for Nf = 6 and 8, and consider a "constant 
Nt" line. As shown in the left panel of Fig. ^ Nt = 6 and 12 lines get to a joint at {g^,N'j) = 
(1.825±0.02, 11.57±0.17), and iV* = 6 and 12 lines at {g{;,N*) = (1.753±0.02, 10.715±0.17), 
indicating the infra-red fixed point with vanishing thermal step scalings. Next we shall 
investigate the critical lattice couplings at Nt = & and 8 for whole range of Nf = — 8. They can 

be well fitted by assuming the functional expression Nf{g-^) = A ■ log B ■ (^gj^ — gj^\]\[^=oj + 1 

giving (c/L^ A^/) = (1-88 ± 0.09, 11.39 ± 0.78) (right panel of Fig. ^. Thus, the thermal step- 
scaling with the use of our lattice critical couplings supports the emergence of conformal window 
near to 12 flavor system, whose (pre-)conformality is now under debate in the recent lattice 
studies 
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Figure 1. (Pseudo) critical values of the lattice coupling g-^ = 10 / f3j^ for theories with 
Nf = 0, 4, 6, 8 and for several values of Nt in the Miransky-Yamawaki phase diagram. Left: 
We have picked up g^ = y^lO/Pj^ for Nf = 6 and 8, and consider the "constant Nt" line 
with Nt = 6, 8, 12. Candidates for the IRFP have been estimated by the crossing points of 
Nt = 12 line with Nf = 6 or Nf = 8 line. Right: The dashed line is a fit with the ansatz 
Nfig^') = A • log [B- {g{: - 5l17V,=o) + l] for Nt = and iV^ = 8 results. 

As indicated by the constant Nt line in Fig. [H the critical coupling is an increasing function 
of Nf for a fixed lattice temporal extension. This behavior is the direct consequence of enhanced 
fermion screening effects due to the larger number of fermion species. We have observed the 
thermal step (or equivalently the asymptotic) scalings for our 13-^, thereby, the enhancement of 
the screening effects gets to a physical significance. The test of the asymptotic scaling has been 
the historical homework since the pioneering finite T study at large Nf by J. Kogut and his 
collaborators [22j and now we have completed it. 

We shall now investigate the Nf dependence of the critical temperature. In order to compare 
a physical quantity such as a critical temperature among different theories with a different Nf, 
it is necessary to introduce a Nf independent reference scale by hand. To set the reference 
scale, ideally speaking, we would like to measure various sizes of Wilson loops at various 



numbers of flavors in Monte-Carlo simulations to obtain the running couplings g in the wide 
range of scales at each number of flavor. In stead of performing such a massive simulation, we 
approximately construct the renormalization flow via the integral of the two-loop beta-function 

Tc _ 1 a~^{9^^) _ 1 f^^' dg 

Wa^) ~ N'tWa^ ~ ^ -9^bo + big^) ■ 

To specify the reference scale M{g-^^^), we utilize our plaquette (tad-pole factor uq) data shown 
in the left panel of Fig. [5J Note that the plaquettes can be regarded as a kind of renormalized 
couplings. Let us consider a constant uq without Nf dependences, for instance uq = 0.9 
in figure, and read-off the corresponding bare lattice couplings at each Nf. The obtained 
gi^{Nf) = ^ylO/ Pi^{Nf ) is used as a reference coupling t/jj®^ in Eq. Q. This procedure imitates 
the scale setting in the potential scheme renormalization, and the use of Nf independent mq is 
motivated by the FRG scale setting method [T^ . 

To be analogous to the FRG study, we should choose uq so as to get a UV M{gi^) 
free from the chiral dynamics. The middle panel of Fig. [5] displays the Nf dependence of 
Tc/M{g-^''^) defined by Eq. P with uq = 0.9. Fitting Tc/Mig-^"^) with the FRG motivated 

ansatz Tc = K\N'f — N f\^~'^^°^^^''^^*f\ we now read-off the lower edge of conformal window 
A'^^ ~ 9.47ib0.02, which is somewhat smaller value comparing to those obtained by the vanishing 

thermal scale settings. In the middle panel of Fig. [21 we find Tc/M{g-^°^) <^ 1, indicating the 
UV nature of the reference scale M{g^^^). 

To get more transparent view for the UV reference scale, we here consider the particular 
reference coupling g-^'^^ - the thermal critical coupling gr^ which makes the reference scale 
M^g^"^ = gr^) be equivalent to in Eq. ([3]): Tc/M{gr^) = 1, giving a typical interaction 
strength at Tc- As shown in the right panel of Fig. [21 the increasing nature of gr^ indicates a 
realization of more strongly interacting non-Abelian plasma at larger Nf as discussed in Ref. [13j. 
The criterion to set the UV reference scale M{g-^^^) at every Nf would be given by the condition 
^^ref ^ gr^i^Nf) for all Nf. We find that uq > 0.84 meets a requirement, while the use of too 
large uq suffers from the strong discritization errors. In practice, we find that the number of 
flavor giving the vanishing Tc/M{g-^'^^) is relatively stable within the range 0.84 < uq < 0.94 
which results in 9.85 > A^) > 9.17. 

4. Summary 

We have investigated the (pre-)conformal dynamics in color SU(3) gauge theories with multi- 
species of fundamental fermions by using the lattice Monte-Carlo simulation. In order to study 
the conformality beyond the fixed number of flavor Nf, we have focused on a reducing chiral 
dynamics at finite T as a function of increasing Nf . We have observed stronger fermion screenings 
resulting from a larger fermion multiplicity at larger Nf. We have investigated a finite T step- 
scaling which follows from the uniqueness of critical temperature (Tc) at each Nf, then the 
vanishing step-scaling signals the conformal dynamics at NJ ~ 10 — 12. Further, motivated by 
the recent FRG based studies [12j, we have examined the Nf dependence of Tc, by introducing 
a UV Nf independent reference scale M{gjJ^^). We have used the thermal critical coupling gr^ 
as a criterion to insure the UV nature of M{g-^'^^) by imposing the condition g^^'^^ <^ gr^{Nf) 
for all Nf. We have found that the number of flavor giving a vanishing Tc/M{g^'^^) is relatively 
stable within the range 0.84 <uq< 0.94, which results in 9.85 > Nf > 9.17. 

As a future perspective, we should measure various sizes of Wilson loops at various numbers 
of flavors, and perform more rigorous scale settings in the potential scheme. It is also mandatory 
to investigate the chiral limit and the thermodynamic limit at large Nf. This, together with a 
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Figure 2. Left: The /3l dependences of the tadpole factor uq at zero temperature with 
the use of 12^ lattice. Specifying a constant uq {e.g. uq = 0.9 in figure), we read off the 
corresponding lattice couplings /3l which are used to define the scale M at each theory with 
Nf. Middle: The Nf dependence of Tc/M where M is the UV scale with uq = 0.9 at each 
theory with Nf. The dashed line represents the fit for data by using the FRG motivated ansatz 
Tc = K\Nj — ]\J'f\^^'^^o/^'^)i'^f) Right: The thermal critical coupling with increasing Nf. 

more extended set of flavor numbers, will allow a quantitative analysis of the critical behavior 
in the vicinity of the conformal IR fixed point. 
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